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Abstract This paper presents a formalism describing the dynamics of a quantum particle in a one- 
dimensional, time-dependent, tilted lattice. The formalism uses the Wannier-Stark states, which are local- 
ized in each site of the lattice, and provides a simple framework allowing fully-analytical developments. 
Analytic solutions describing the particle motion are explicit derived, and the resulting dynamics is studied. 
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1 Introduction 

A recent revival of interest on the old problem of elec- 
tron motion in perfect, non-dissipative, periodic lattices 
is due to the advent of laser cooling techniques, that 
have opened the possibility of realizing the interaction of 
atoms with perfect, non-dissipative, optical lattices. Op- 
tical lattices are a consequence of the displacement of 
atomic levels resulting from the interaction with light 
0]. A far-detuned standing wave formed by two counter- 
propagating laser beams (of wave- vector k^) is perceived 
by an atom as a one-dimensional potential V , whose strength 
varies sinusoidally in the space: V oc cos{2kLx) (where x 
is the spatial variable). Standing waves are the model po- 
tential considered in the present work. The main source 
of dissipation in such lattices is spontaneous emission, 
which can be reduced in a controlled way by increasing the 
laser-atom detuning. By scanning the frequency of one of 
the beams forming the standing wave with respect to the 
other, one can generate an accelerated potential. In the 
(non-inertial) frame in which the standing wave is at rest, 
an inertial constant force appears, producing a "tilted" lat- 
tice. This technique can also be used to modulate (tem- 
porally) the position or the slope of the potential. Tilted 
optical lattices generated in this way have been used for re- 
cent experimental observations of Bloch oscillations both 
with atoms (Sj or Bose-Einstein condensate jBj. Wannier- 
Stark ladders |Zj and collective tunneling effects [H] have 
also been studied with such a system. 

A very convenient theoretical framework for studying 
atom dynamics in a tilted lattice is that formed by the 
so-called Wannier-Stark (WS) states which are the 
eigenstates of a tilted lattice |l()llllj . In particular, we 
have studied in a previous work the quantum dynamics 



in a time-modulated, tilted lattice using such basis |12l 
111) . In the present paper, we present exact solutions of 
the equations of motion obtained in Ref. ^I] that describe 
the atomic center of mass motion. Sec. |21 briefly reviews 
the main properties of WS states leading to a simple de- 
scription of Bloch oscillations; Sec. 13 studies the coherent 
dynamics in a modulated lattice, leading to a general equa- 
tion of motion; we derive in Sec. 0] an exact solution for 
the equation of motion which is studied and characterized 
in Sec. [SI 



2 The Wannier-Stark basis 

Let us first introduce natural units in which lengths are 
measured in units of the lattice step d {— Ai,/2, Al = 
27r/fci being the laser wavelength), energy is measured in 
units of the "recoil energy" Er = h^k\/ {2M) where M is 
the atom mass, and time is measured in units of h/ En. 
The Hamiltonian corresponding to a tilted lattice is then: 



p2 

Ho = h Vb cos(27ra;) + Fx 

2m 



(1) 
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where m = 7r^/2 is a reduced mass, F is a constant force 
measured in units of Eji/d, the momentum operator in 
real space is p = —i{d/dx), and h — 1. 

As shown in Fig. ^ the symmetries of the tilted lattice 
suggest that the eigenenergies shall form "ladder" struc- 
tures separated by the "Bloch frequency" lub = F (or 
ujb = Fd/h in usual units), the so-called "Wannier-Stark 
ladders" [9 . Each ladder corresponds to a family (labeled 
m) of eigenfunctions ipnm(x) (the Wannier-Stark states) 
centered at the well n, and thus spatially separated by an 
integer multiple of lattice steps. Inside each family, WS 
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Figure 1. The tilted lattice and Wannier-Stark states. The 
WS state labeled (a) is a delocalized state of the continuum. 
States labeled (b) and (c) are lowest-energy states localized, 
resp., on the sites n — —1 and n = 1 and are identical except 
by a translation of two lattice steps. The energy levels drawn 
at the right of the plot show the Wannier-Stark ladder. 



states are invariant under a translation by an integer num- 
ber n of lattice steps, provided the associated energy is also 
shifted by the same number of wb, i.e. 



and 



Enm — Eom + flUJ B 



(2) 



We consider a spatially limited lattice, extending over 
many periods of the potential, enclosed in a large bound- 
ing box ^. Although the presence of the bounding box 
breaks, strictly speaking, the above symmetries, it changes 
only very slightly the properties of "bulk" states, and the 
eigenenergies and eigenstates obtained numerically dis- 
play, as long as we stay far from the bounding box, these 
symmetry properties to a very good precision. The exis- 
tence of WS states has been evidenced experimentally in 
1988 in a semi-conductor superlattice |13II14| . and 1996 
with cold atoms in an optical lattice [7]. 

We shall consider here parameters such that WS states 
are strongly localized inside a given lattice well, and that 
the atom dynamics can be described to a good accuracy 
by the lowest-energy WS state in each well (as the states 
noted (b) or (c) in Fig. ^. We thus drop the family 
index m, and note the WS states by (/3„(a;), with associated 
energy En = ucub (setting Eq = 0). 

The system dynamics is described by projecting the 
atomic wave function on the WS basis: 



(3) 



^ In the absence of a bounding box, Wannier-Stark states are 
metastable states (resonances). 

^ This is equivalent to neglect Landau-Zener inter-band tran- 
sitions in the usual Bloch-function approach. 



with Cn{t) — c„(0)e^*"'^^*. Dynamical observables can 
be easily calculated. For instance, the mean value of the 
atomic position operator, {x)t — {'P \X\ ^) is: 

{x)t ^^Xnn \Cn{t)\^+Y^ (X„™< (t)c™ (t) -f C.C.) (4) 



where Xnm = {^Pn \X \ ipm) are constants depending only 
on the parameters Vq and F of the tilted lattice. As we 
are considering WS states localized in the wells, we can 
keep only nearest-neighbors contributions (|n — to| < 1) 
and derive a simplified expression: 



{x)t 



] Cn(0)*c„+i(0)e 



+ C.C. 



(5) 



where x — ^ Xnn |cn(0)|^ is the mean position of the wave 
packet and Xi = X„(^„+i){^ ^oi = -'^o-i = ^n(n-i)) is 
independent of n ^. Eq. JSJ describes Bloch oscillations, 
in a simpler and more intuitive way than the usual Bloch- 
function approach |[^15| . The amplitude of the Bloch os- 
cillation is proportional to Xi , and grows with the overlap 
between neighbors WS states, i.e. it increases as the slope 
F and lattice depth Vq decrease. The physical origin of the 
Bloch oscillations appears very clearly as an interference 
effect between neighbor sites, since c* c„+i is the quantum 
coherence between the sites n and n + 1 Bloch fre- 
quency is seen to be just the fundamental Bohr frequency 
of the system. 



3 Coherent dynamics in a modulated lattice 

We now investigate the case in which the system is "forced" 
with a frequency close to its natural frequency lob- We 
therefore add a sinusoidal component of frequency uj [i.e., 
a term of the form Eq sm{ojt)x] to the Hamiltonian Eq. QJ, 
with UJ = ujb + S and \6\ <C co^cob- The WS basis again 
leads to an analytical approach TT. We limit ourselves 
to modulations that are smooth enough in order to avoid 
transitions to delocalized or to localized excited WS states, 
ensuring that the decomposition over the lowest-energy 
state of each well [Eq. (0)] remains valid for all times. 

The coefficients c„(t) can then be obtained by report- 
ing Eq. Q into the Schrodinger equation 



. d<F{x,t) 



p2 

2 m 



+ Vo cos(27ra;) -t- Fx — Fqx sm{Lut) 



(6) 

which produces the following set of coupled differential 
equations for the c„(t): 



i{t) = —inujBCn{t) + iEq sn\{ijjt) 



X„ 



By virtue of the translational properties of the WS 
states, Xnn+p (p 7^ 0) does not depend on n: Xnn+p ~ 
J tp1^{x)xipn+p{x)dx = J" iPq{x — n)xipp{x — n)dx — 
J LP(){x'){x' + n)ipp{x')dx' — Xop = Xp, where we used the 
orthogonality of the WS states. The same kind of calculation 
leads to X„n — Xoo + n. 
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where c„ = dcn/dt. Neglecting temporarily the couphng 
between different WS states, (i.e putting Xnm — Xnn^mn), 
the amphtudes are obtained simply as c„ — exp[i0„(i)], 
with the time-dependent phase 



Equation (|12|l is similar to a "dipole coupling" between 
sites n and n ± 1, where J?i plays the role of a Rabi fre- 
quency ^ . In the next section, we will find an exact solution 
for the above equation. 



(t)n{t) = -nujBt - ( — ) (-'^oo + n) cos{ujt) (7) 



where we used the identity X„„ = Xqo + Writing 

c„(i) = d„We^^"'*\ (8) 
the new amplitudes dn are governed by: 

dn = iFo ^ Xnmdmit) cxp {i [(pm{t) - 0„ (i)] } sin(tjt ) . 

(9) 

After Eq. Q, the phase difference (l>m{t) — <f>n{t) is: 



(l)m{t) - 0„(t) = [n - to) 



UJBt 



Fo 



cos{ujt) 



where we used Eq. ((2Jl. Eq. @ can be recast as: 



dn — iFo ^ Xpdn+p ^ 



g ipujBt^—ip{Fa/uj)cos(Lot) 



sm(ujt) 



p^o e ^ 

X [exp {i [{£ + l)uj - pllib] t} 
— exp {i [{£ — 1)lu — pujs] t}] 



(10) 



where Xn 
tion of t 
formula: 



X„ 



(p 7^ 0)121 Jn{x) is the Bessel func- 



^p — ^^n{n+p) 

tion of the first kind, and we have used the well-known 



-f-oo 



-iz cos[Ljt) 



(11) 



The coupling parameters Xp rapidly shrinks to zero as 
p increases and we can keep, to a good accuracy, only 
the contribution of the neighbor site (p = 1) (for in- 
stance, Xi = bx 10^2 ^ ^2 = 8 X 10"'' for Vo = 4.5 
and F = 0.5). Moreover, the sum over the harmonics of 
w (i.e. t) is also limited to a few terms close to £ = 
[typically, irnax ~ pFo/uJ ~ 0(1)] Finally, we keep only 
the slowly varying terms in Eq. (|10|) : since 5 — lo — ujb 
is assumed small {\5\ <ti lu,lub), on can retain, to a very 
good accuracy, only the terms which oscillate as exp{±iSt) 
(the fast oscillations give a vanishing contribution on the 
average). To the leading order, this gives: 



dn{t) = Hi [dn+ie'^* ~ dn-ie~ 



■iStl 



(12) 



where 



Hi = 



FqXi 



UJ J \ UJ 



^OjXiJi ( — 



(13) 



The Bessel function value Jn{x) becomes small for |a::j 



4 Exact solution 

It is interesting to consider the complex amplitudes d(k, t) 
defined in the reciprocal space: 



d{k,t)= J2 ^«We"' 



n— — OQ 



dn{t) = — / d{k,t)e-''"'dk. 



(14) 
(15) 



In these expressions, dn{t) are the Fourier coefficients of 
the continuous and periodic function (period 27r) d{k,t), 
and k plays the role of a quasi-momentum. Eq. p2|l can 
easily be translated to the fc-space: 



d{k, t) = 2if2i sin((5t - k)d{k, t). 



(16) 



The quasi-momentum distribution is obtained after a straight- 
forward integration: 



f 21^1 

d{k, t) = d(k, 0) exp < — - — [1 — cos((5t)] cos(fc) 

r 2z/?i . . 
X exp < — - — sm(ot) sm(— fc) 



(17) 



and, using Eq. (|ll|l and the equivalent relation e'^ ^''"^ 



-iik . 



d{k, t) = d(k, 0) E Jp ( ^ [1 - cos(Jt)] j 

pX ^ ^ 



(18) 



Setting £ = p + q: 



d{k, t) = d{k, 0) E Jp ( ^ [1 - cos(5i)] j 

p,q ^ ' 



(19) 



and using the Bessel addition theorem (see Appendix 
we transform the above equation into: 



d(fc, t) = d(fc, 0) E Jg(Q)e"*«''e*' 



qSt/2 



^ Eq. predicts that the coupling vanishes if Fq/uj co- 

incides with a zero of J\{x). However, the first root of this 
function arises for x ~ 3.8, that is, for Fq ~ 3.8a; = 3.8-^, 
which cannot be considered as a smooth modulation, falling 
outside the range of validity of the approximations leading to 
Eq. (O. 
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Figure 2. The gray-level plots show the temporal evolution 
of the probability of presence of the atom, initially localized 
in the site m = 0. Plot (a) shows the resonant behavior: the 
wave packet spreads linearly with time. Plot (b) displays the 
non-resonant "breathing" behavior. {S is indicated in units of 
4f2i). 



where 



g(t) = ^sin(-). 



It is now straightforward to obtain the coefRcients d„(t) 
by applying Eq. (|15|l : 

d„(i) = ^rf„+,(0)J,(Q)e^«^*/2 (20) 
= ^d,„(0)J™_„(g)e'("-")**/2 

m 

This remarkably simple result shows that amplitude of 
probability of finding the particle in the state n at time t 
is simply the sum of the contributions of the initial pop- 
ulated states m, c?m(0), with a time dependent weight 
Jm-n (Q) e*^™"")'^*/^, that depends on the distance n — m 
between the sites. The presence of the evolving phase term 
gi(m-n)<5t/2 evidences the quantum-coherent nature of the 
dynamics. 

The most relevant contributions to the above summa- 
tion are due to sites such that jn — m| < |Q(t)|. Roughly 
speaking, the dynamics of site n is correlated to the sites 
situated inside a range An « |(5(i)| around n, which, out 
for resonance, is roughly |4i7i/(5|. Close to the resonance, 
this range is An « 2 |i7i| and the dynamics on site n 
will "mix" contributions of larger and larger numbers of 
sites as time increases. To give a simple example, consider 
the case d„i{0) — Sm.o- One then finds {6 ^ 0): 

\dn{t)\^ = [J„(g)]' 

which shows that the n^^ site is "filled" and "emptied" 
in a periodic way with a period 27r/(5. This "breathing" 
behavior is shown in Fig El 



5 Wave packet dynamics 

The dynamics of a wave packet can be characterized by 
calculating mean values of observables. This can be done 



analytically in the present framework. One can easily de- 
duce from Eq. (IH) the mean position: 



{x)t = ^X„„ \d„it)f+XiJ2(dn{trdn+i{t)e 



i<l>(t) 



+ c.c. 



(21) 



where 4>{t) = 4>i{t) — 4>o{t) = —ujBt—{F/uj)cos{ujt). Using 
the relation X„„ = Xqq + n = Xq + n, a, straightforward 
although somewhat lengthy calculation (see Appendix Ib)| 
produces: 



{x)t = {x)t=Q + 



Q{t) 



-i&t/2 



i4>(t) _ gi0(o) 



X Y,di{Q)dn+m + c.c\ 



(22) 



The mean position will evolve with time only if the ini- 
tial packet displays site-to-site quantum coherence, i.e. if 
d* (0)(i„+i(0) ^ (this is also a condition for observ- 
ing Bloch oscillations, as pointed out in SecEJ. The evo- 
lution in time may display both a fast oscillation at the 
frequencies lub and m (and its harmonics) and of a slow 
oscillation at the frequency 5. Note that the weight ra- 
tio between the slow and the fast component is roughly 
fiil{5Xi): the slow oscillation is dominant close to reso- 
nance. 

The meaning of the above result can be more clearly 
appreciate by considering the simple case in which the 
initial wave packet spreads over N sites and presents a 
constant phase difference e^^° from site to site: 



d„(0) 



N 



(23) 



Substitution into Eq. H22|l leads to: 

{x)t = {x)t=o + {~Q{t) cos(fco - St/2) 

+ 2Xi cos [fco — o^st — (Fq/lj) cos(cjt)] 
- 2XiCOs(/co - Fo/u))} 



(24) 



If (5 = 0: 



{x)t=o — — 2i7i cos{ko)t 

+ 2Xi (cos [fco — LUst - 
- cos (fco - Fo/llib)) = 



{Fo/ujB)cos{ujBt)] 
vgt + fit) (25) 



where Vg = — 2i7i cos(fco) is the group-velocity of the wave 
packet, and f{t) a purely oscillating function. This result 
shows that the main features of the resonant dynamics are 
determined by fco. If fco = ±7r/2 (phase-quadrature from 
site to site), Vg = and there is no global motion (but 
there is an oscillation and a spreading of the wave packet). 



If fco 



vr, Va 



2J7i and the global motion is a climb up 



along the slope of the potential with the maximum speed 
2i7i. In the case fco = 0, where Vg = — 2!7i the atom 
climbs down the slope of the potential with a constant 
maximum speed: there is coherent transfer of energy from 
the modulation to atom, thanks to the particular phase 
relation between neighbor sites. Contrary to the motion 
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of a classical particle, the speed \vg\ is independent of the 
sense of the motion: the wave packet climbs the slope up 
or down at the same speed. 

Let us finally note that the spreading of the wave 
packet {x'^)^ can also be obtained analytically. However, 
this leads to heavier formulas, that the interested reader 
can find in the Appendix IbI 

Before closing this section, let us note that other res- 
onant behaviors can be observed for uj = qtOB (<? integer). 
From the general expression of Eq. one finds for a 

coupling between sites of the form (n — > n± q), for which 
the same techniques can be applied. For example, resonant 
interaction with uj = 2lub, leads to: 



dn{t) = J72 [d„ 



dn-2] 



(26) 



with J72 = wA'2 Ji(-fo/^)- Numerical calculations show 
that X2 <^ Xi. Around this resonance, interesting kinds 
of coherent dynamics can be observed 



6 conclusion 

We have developed an analytical approach to the dynam- 
ics of a wave packet in static and time-modulated tilted 
lattices, and provided exact solutions describing the atomic 
motion. It is worth noting that the present description is, 
in its principle, independent of the shape of the lattice, 
provided it presents localized states. It is therefore gener- 
alizable to non-sinusoidal lattices. 

The major experimental difficulty for observing the ef- 
fects described in this paper is the creation of the initial 
atomic coherence. There are various solutions for that. 
One is to cool the atoms to sub-recoil temperatures, so 
that their de Broglie wavelength is of the order of a few 
lattice steps (5-10 lattice steps is typical) [5]; the light 
potential is then turned on adiabatically. Another way of 
creating spatial coherence is to start with a Bose-Einstein 
condensate, whose spatial coherence length is up to hun- 
dreds lattice sites (HHEI- Note that in the latter case the 
system obeys a nonlinear Schrodinger equation, but the 
present approach can still be generalized to such case T^ . 

The detection of the coherent dynamics described above 
is simpler in the momentum space, using velocity-sensitive 
Raman stimulated transitions jElElEl- The coherent mo- 
tion of atoms, e.g. when they climb the slope of the mod- 
ulated potential (Sec.|2I), correspond to a speed of about a 
recoil velocity, i.e 3 mm/s for cesium, which can be easily 
detectable by Raman stimulated spectroscopy. It would 
be nice to observe the coherent motion also in the real 
space. This seems to be very difficult, because the spatial 
amplitude of the motion is very small, of the order of a 
few lattice steps, as compared to the atomic cloud which 
extends over hundreds of lattice wells (for atoms cooled in 
a Magneto-Optical trap). The coherent dynamics is thus 
detectable in momentum space, but very hard to see in the 
real space, at least under usual experimental conditions. 
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A Bessel addition theorem 

The addition theorem for Bessel functions of integer order 
states that: 



^ Jp(ro)e-^^'°Jp+,(ri)e^(^+''"'i = MQ)e 



iqO 



(27) 



where ^ = ~r*i 



r 



R e''-', and r? = rie" 



ro = roe 



iBo 



For Eq. (O, one sets m = [1 - cos{5t)] , Oq = -7r/2 ri = 
sin(5t), 6*1 = 0, and finds quite straightforwardly 



and 



4^?l . St 



Eq. lion takes the following form: 

d{k, t) = d{k, 0) Jp [1 - cos{5t)]^ 



p. 9 
2Qi 



ip^ -iqk 



= d(fc,0)^ J,(Q)e~''"=e'«* 



B Mean values 

Consider first the first term on the R.H.S of Eq. (I21t : 
{x)['^ = ^(Xo + n) ^ J„.(Q)d;+„(0)e-™*'/^ 

n m 

X ^ J,(Q)d„+,(0)e"**/^ (28) 
e 

Changing variables so that q — m~£,p — n + £ leads to: 



{x}r' = ^ {Xo+n)dn+e{0)d*,_f.i+q{0)JiJe+ge' 



iqSt/2 



= ^ {Xo+p-£)dp{0)d;+^{0)JiJe+,e-"'''^'' 

where J„ = J„ (Q), Q = ^ sin(f ). Eg . I?7l implies here: 

J2 MQ)Ji+q{Q) = JM = Sq,o. (29) 

e 

and using also the Bessel recurrence: 



(Ji+iiQ) + Ji-i{Q)) 



(30) 
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one gets then 



(|e-*'/^X^d,(0)d;+i(0) + c.c. 



sub It is simpler to show, by the same kind of reasoning, that 
the second term in Eq. 12111 is 

{x)f^ ^XiJ2 (dp(t)*dp+i(t)e'^(" + C.C.) 
p 

= Xi ^ (dj,(0)*dp+i(0)e'^(" + c.c.) , 
p 

so that the mean is: 



^dp(0)*dp+i(0)j +c.c.| 



with 



p 

+ Xie'^(°) Yl (rfp(0)*rfp+i(0) + c.c.) 
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One can calculate, by the same kind of technique, the spread- 
ing of the wave packet: 



Xo+p+- ]e- 



iSt/2 



(2) ( mt) MO) 



d;(0)dp+i(0) +c.c.| 



+ X S [(d;(0)dp+2(0)e-'^* + |d,(0)r + c.c.)](31) 

p 

where x'^^ = J ipQ{x)x^tp\{x), and 

V 

+ Xf e"^(«) (rfp(0)*dp+i(0) +C.C.) 



with = J <^S(2;)a; H>o[x) ~ J 'Pnix)x <Pn{x) — 2nXo — n . 
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